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Magneto-hydrodynamic Oscillations of a Conducting Liquid 
Mass Rotating in a Uniform Magnetic Field 


By Tomikazu NAMIKAWA 
Geophysical Institute, Kyoto University 


(Read Oct. 28, 1955; Received May 25, 1956) 
Abstract 


Considering Magneto-hydrodynamic oscillations of an infinitely long cylin- 
der of electrically conducting liquid medium rotating in a uniform magnetic 
field, we see that the periods of oscillations are much different from that of 
non-rotating cylinder, when Coriolis force is much larger than the electro-magnetic 


force. 
We obtain two periods 7,, 7, for one mode of oscillation 
T\>TM>T: 


where 7, is the period of non-rotating cylinder. 
1. Introduction 


The oscillation of non-rotating conducting liquid under the influence of the 
magnetic field has been worked by many authors [1], [2], [3]. In this paper we study 
the same problem taking rotation in consideration. The Coriolis force, which usually 
plays no significant role in physical phenomena of laboratory scale, is found not to be 
neglected in phenomena of cosmical scale. 


2. Fundamental Equations 
Consider a homogeneous mass of infinite electrical conductivity partaking in 
rotation with a constant angular velocity @; and let there be present also a uniform 


magnetic field of intensity Hy applied from outside. The fundamental equations can 
be written as 


curl A=4n7 : curle = — (1) 
e+vx(M\+h)=0 ; H=H+h ; (2) 
div h=div v=0 (3) 
dv 
and oof oy. + (O+P)D-+-28x vo} =ix(Ho+h)—PP (4) 


where i, e, h, 0, P and v denote respectively the electric current density, electric field, 
magnetic field, density, pressure and velocity. It is also assumed that the magnetic 


permeability is unity, and that there is no non-electro-magnetic force which depends 
on time. 


An elimination of @ and e gives 


( 9 ) 
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ve = (Ay-7)v+curl(v xh) (5) 
dv 1 1 
a, = ~~ (Ao) h4+20x2—pd——— Axcurlh+vxcurlv (6) 
ot 40 Arp 
and re Pe 1 (Hyrh)+ © (7) 
09 4200 


In the following we take only small oscillations in account, so that the seconc 
order terms may be neglected. The equatiohs (5), (6) and (7) take the form 


te =(H):7)v (8) 

= Pia (Ho: 7)h+2v0xQ—r¢ (9) 

ot 400 : 

P =1{ P+ : (Hy) (10) 
Po 4n 


Taking the curl of equations (8) and (9) and making use of solenoidal property 
of v, we get 
0b 


= : a 
S: (H)-7)@ (11) 
Ow 1 
: A ‘ 3 12 
a PIR 7)b+2(2-7)v (12) 
where b=curlh and w=curlv (13) 


From (11) and (12) we obtain 


o? 1 Ov 
“py = : 14 
a2 A Pla P) lo 2(2-7) at : ) 


Remembering that, in view of the solenoidal cearacter of v and h, 
curla=—p’v and curlb=—pj’*h (15) 
the operator 7? acting on a vector, being understand as 
y?=7V-—-V xVx ; 


we take again the curl of equation (12) and obtain 


= 0 7°v=2(2-7)w—-——(Hy-p ph (16) 
Ot 4709 
Differentiating this equation (16) with respect to time and making use of equation (8), 
we get 
0? Ow 
H,-7)? |p?v= —2(2-7) —— (17) 
[on Fagg (Her? |p'0= 200-7) 
From equations (14) and (17) we find 
+ *17°v=— : 18 
loan ~ 44; (He 9' | 0 = 4(Q-p)o (18) 


Using relation (3), a divergence operotion on equation (9) gives 


y°¢=2 div(vx @) (19) 
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Measuring the length in unit of a, typical length of the size of the conducting 


mass, equations (8), (18) and (19) take the form 


Oh (MM. 
Ot a a rp 
o2 u 27 7 ,0"0 
— —~*e = =—4 Q- lid . 
ee (= r) [rw (Q-7) a 
p’°¢=2adiv(vxQ) 
wee 
where Or Peas 


Equations (20), (21), (22) and (23) are the basic equations of this theory. 


3. The Case when H, and @ act in the same direction 


(20) 


(21) 
(22) 
(23) 


When the axis of rotation and the direction of the impressed magnetic field 


coincide, we choose a set of cylindrical coordinate (p, g, 2) in which the z axis lies 


along the direction of Hj and @. 
We can separate the variables by assuming 
D« ettzeime pire 


then equations (20), (21) and (22) have the forms 


h= Hy ! D 
ah; 
p’v=—K*v 
V°¢=2 a2 div(vxk;) 
or p9=209 1] (ove) — S| 
eL O% dg 


where K=—ap— and k, is the unit vector in the z direction. 
en A 


4. Solution of Equation (26) 


Let the scalar function,” be a solution of the equation 


Pp +K*p=0 


we now construct three independent vector solutions of (26) as follows [4] 
L=rd, M=rxko, N= 7xM 


Since k, is a constant vector, it is clear that M can be written also as 


M=Lxky= (7x 


(24) 


(25) 


(26) 
(27) 


(27’) 


(28) 


(29) 


(30) 


The vector functions Z, M and N have certain notable properties that follow directly 


from their definitions. Thus 


VxXL=0 p-L=p*p=—K%” 
Whereas M and WN are solenoidal 


y-M=0, y-N=0 


(31) 


(32) 
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As velocity and magnetic vector are purely solenoidal, the expantion is made in 
terms M and WN alone. 
1. Axial symmetric solutions 


19 = Jol /K*—Ep) ee! ie os) 
il Oj, 
=A 02 0 pilzpiat 
One ig 
Vg= 2B 6 2 elteeir= (34) 
Bn 
v, =A, = S Juatl ‘on! 
2. Asymmetric solutions 
9)" =J..{ / K?—Po) ei * ett ett (35) 
e, =»| 4" al A+B, mimi. |etmretite'™ 
dm K 0 0 
Ve = —>} [amie Vy. bE Pies = etm? pilzpinrt (36) 
im K op Op 


Ke 


‘es sacar . Oe 
J,e' ™ P pilzpi Xe 


v,>= DA re 
dm 


5. The Magnetic Field Associated with the Oscillations 


For points within the conducting mass, the magnetic field is given by (25), (34) 
and (36) as follows, 
1. Symmeritic solutions 


_ Ah PY haw al? ake ilzpidt 
hy at YA ap” 
= ig of OJ wlzpidt (37) 
~— a UB A 0p ee 
_ of (K?—P) ilzpirt 
mee a BA, Ki rene 
2. Asymmetric solutions 
_ i OJ, meml a im® ptlzpidt 
ash) Sicacaee bie iondibs dase © 
_—-Mh nem des m L a im@ ptlzptrt (38) 
hy = hal a | a i he ee 
h.= Ho yy 4 ml KEKE) pime gitegiat 
= a im Ki 


6. The Variation of Pressure in the Mass 


Considering the resolute of equation (9) in the z direction we find that ¢ is of 


the form 
Symmetrical solutions 


K?-P 
¢°=2a2>3A,° kK? YF aed deed (39) 
Z 


Asymmetrical solutions 


g=2 a2rA." ne BI ewe a (40) 


Magneto-hydredynamic Osawhatietil of « Conducting Liquid Mass wine 


gi sbsm ef Mbiaasoxes att fobioholee udinuad a8 yWOIH Nie aot e h si nots 
4 _ Mesa, equations (8), (18) and (1M take the form onols Vi tos 
j ah (.\, en me eh 
«) OF oy A—RY Wo — m5 
re /g *Fa. a tid 
| ae ~ Sn} fal “WRK, 
(be) . r“¢ =e Gere *art— = eB <G (22) : ; 
where a= a ¢ > (2) 
v tageigy = = 2 Renae 7 ag 
Equations (20), (21), (22) and (23) are the Basic equations of of thie eee oA 
(2) 2. The Com, when fF got 4 im the same direction « 
When the @: rai “age ee aici, rene pM pressed magnetic field 
coincide, we chhose a set of lofindrical BbrPrate S¢,.2) in which the # axis ies 
© AORe the directiba OOM wHs Bag yh x Ag <a 


an mae oi 
/ 
f 
yt 
7 : 
i 


We can sebarate the vartables by assuming, > : o 
3Ab = — 
a : Me ea a Ce ’ 
Rs then ©qosiiibindMbiat? a0 dai bbtenGdKA SIGE sitemeaM od? .2 a 


ae 


8) (28) ed vig a boa siting Sie 22em yaitsubaoo ai aidsiw 2tniog 10 ee. ox 


. Figs = iy emoltis a tie ee f 
ea eer tr ee 
a P'¢= Pale kT ae ney | er) << 
a att — 
where a / 


= 


: - dig etng | 

- Kam - nd &, is the 0 het Rhe z direction; 

i a? chp cae age. aA 

* 4. Solution of Equation nln Pr Pron Gi 

; Let the exalar fu 6 ae 

z ms Way “ S. REE ee igual es 
. 4 bey + and ee 


‘ BB ow constridt meets Totbut *haaadaihersan {4} — 
be dit en et ee 
a an mer lagea en : 
es ® tans bah sw —— to-stuloesy only anata? 


a Phe veetor functions Z, Mf anc id A’ have certain notable | nant tig 
fom their definitions Tins ge 
(RE) 1 me” 


= Pxbod ar x 
- Whereas E and W are soleno ida} 
Ob) a” 


. ‘eo A, a oa 


7% 


ae ee a ee ' 
ie el a 


“St See Beg CTS) eitsups mont to 
We peer eek.) angituloz lsorgammye 


reer pres: a Rn ae ee ee ae 
2noduloe [poittorinnyeA 


. L on re «a paca 
a aM Ateous. =a 
os Bie the Rapti — 9 tug sw enoiisups wot svods mord 

4 . | "Ds ae wa ’ —< : 
a “K-=".6 “Bns 
= geeA Inova oirtoolH oT 5 a 

oa “Thad ay > "Ri = 
= 198 ow Ww) noiEups mort 
a. 

=, — 
iii. atl os 9eU aon 


Le “ ti SS af 


where ani: lime | <p 


e 


Magneto-hydrodynamic Oscillations of a Conducting Liquid Mass Rotating 101 


or from equation (27’), we find 
Symmetrical solutions 


o° = —2a2>SB,° gees lee (41) 
Asymmetrical solutions 

o=—2a99B" © oP Tueltemrein | (42) 

ee kK: 
From above four equations, we get 
os 0 

or Fea } (43) 

and B= —A,™ 


7. The Electric Current Associated with the Oscillations 


From equation (1), we get 


i=— curth (44) 
4na 


Making use of the relations (28) and (29), the electric current density are written as 
follows: 


Symmetrical solutions 
Hy _ B,° wl? Of, ei tzeirt 
dp 


i= = a Z 
4na* 1 A 
ig= stiles gated epban (45) 
4na* i a 
tye ee a2z pide 
a 4ra* SB, O he 7 
Asymmetrical solutions 
- __ WP OJ, tml Bi Jo | sme stmeetht 
tp= re =| ae oy Op +B," ma =e ée e 
H, P-m OJ 1K OJ 
= A m m m a ™ |oim®? pilz pirt 
, 4na? i. =| as 0p te A Oo i ve cm 
en Hy Ae —P) im? piles pidt 
1, = “Gra? et ya paige g 


8. The Electric Field Associated with the Oscillations 
The electric field within the mass is given by 


e=—vuxH. 
Symmetrical solutions : Hs? 


ep =H,>)B,° Jp etttetrt 
? Op 
‘ lO (48) 
Co = Hi 4,0 6. CL0 pitz gine 
° oA, K he e 


¢.=0 
Asymmetrical solutions 


mh lm J, 0) 
ép=AM)>}| A,” "4 Bm J™ | eimevitzeine 
e oe | K é ap Je e e 


il OJ, 


: ' 49) 
"41 Bow was im? p tle pi rt ( 
K dp + 2 tm - k ene 


eo = Hox »[ a 


é,=0 
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9. Infinitely Long Cylinder 


We picture to ourselves an infinitely long cylinder of uniform circular cross-sec- 
tion of radius a, rotating with constant angular velocity @ about the z axis, along 
which a uniform magnetic field of intensity H) acts. 

1. The period of the normal modes of vibrations. 
a) Rigid bounding surface. 
When the bounding surface is rigid, the bouhdary condition is 
ters at” p=1 (50) 
z»p=0 
1) Symmetrical solutions 
In this case as the g¢ component of the external magnetic field is zero, we have 


so, the period equation for the normal mode of vibration is 


Oy =F) =o or J,(,/K?—2)=0 (52) 
fo) 
where _ 201 (53) 
~ we 
a. 
a 


When a; represents the roots of this equation, we have 


K?—-P=a; (54) 
Making use of relation (53), we get 
Q ul? 
24+2—A——.=0 (55) 
Art 1, 2 
where bs: 
2 
L= oF + (56) 
Solving equation (55) we find 
A= +A 4/1+R?+R) (57) 
where P A 
ples Reet (58) 
dot s* z Liu Udo 


R corresponds to the ratio between the Coriolis force and the electromagnetic force. 


From the relation 


ax 2 4 (59) 
<i ea d 
we have 
T=Tp(./1+R? +R) (60) 
eee 2na 
where T is the period of oscillation and N=, (61) 


repsesents the period of oscillation, when the cylinder is not rotating. 
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From the relation (58) we may write 


7d in jos (62) 
| do | Tt 


which indicate that Coriolis force is significant for a vibration with a period T com- 


parable to the period of revolution rt, of the cylinder. 
If R>1 (60) becomes 


ey 
rik 63) 
aT ToR{(1+ 9 x)=" ( 

We have 

T,=2RT 6 (64) 
and 
1 

~ 65 
T; 2R Ts ( ) 


2) Assymmetrical solutions 
In this case the period equation for the normal mode of vibration is 


1. Ge mIn=0 (66) 
As 
Yu P)_ KR Ff ee e: zi} (67) 
we have 
1/KP—PJy_s(./K°—P)—m (K+ DJn(/K?—P)=0 (68) 
Putting 
JK?-P=x (69) 
we get 


K?=24+2 or K=J/x+P 
Equation (68) becomes 
IxJn\(x)—m(L4J/x*®+P)],,(x)=0 (70) 


In case when m=1, equation (70) takes the form 


lxJo(x)—(L4J/x*+P)],(x)=0 (71) 


b) free bounding surface adjoing vacuum 
At the boundary of the cylinder 
D¢ _ = 
Di =0 at p=1 (72) 
and to the first order this equation implies that 


¢=0 73 
So we have from (39) and (40) (73) 


Ini K?—-P)=0 (m=0, 1, 2, 3,-+--) (74) 


But other boundary conditions requires 
Symmetrical solutions 
AEE Sits P) _9 


(75) 


a= wre 
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Asymmetrical solutions 


la CF 
gets = as Pa T /K*—P) =0 (76) 


There are no roots common to J “=0 and J,,=0. 
0 


We have no solution in a cylinder having free bounding surface adjoing vacuum. 
10. Application to the problem in the earth’s core 


Taking /=2z (corresponding to the wave length of core’s radius), a=3x108cm, Hy= 
3 gauB, o)>=10 gr/cm®, @=7.27x 10-5, @,=3:83 
we have 

U=0.27cm/sec, 1,=1.17, R1.1X104, Ty=11.1x 108sec. 


Making use of these values we obtain from (64) and (65) the value of periods 


T, 7.7 X 10° years 
T.~14 hours. 


These values may be taken for representing the order of magnitude. 
11. Conclusion 


The discussion of the previous section has shown that the Coriolis force plays 
an important role for magneto-hydrodynamic oscillation in the earth’s core. 

As the model is very rude, the results are not directly applicable to the real 
case, but it may be worth while to notice that in a rotating mass, a mode of large 
scale motions has two periods, one is longer, and the other is shorter, than that of 
non rotating case. The writer will report on the oscillations in a rotating sphere in 
the near future. , 

In concluding this paper, the author wishes to express his thanks to Prof. M. 
Hasegawa for his constant interest in and encouragements and to Dr. M. Hirono and 
Dr. K. Nagashima for their valuable discussions. 
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Anomalous Increase in Thermo-Remanent Magnetization 
of Ferromagnetic Minerals.* 
— Magnetic Interaction between Different Constituents 
in Ferromagnetic Minerals — . 


By T. NAGATA and M. OZIMA 
Geophysical Institute, Tokyo University 


(Read May 13, 1956: Received May 31, 1956) 


Abstract 


Magnetic interaction between different constituents in each grains of ferro- 
magnetic minerals in rocks is studied with the aid of X-ray, microscopic and 
chemical analyses as well as magnetic measurements. Since by means of suitable 
heat-treatments the degree of intergrowth of a constituent into others in mineral 
grains is changed markedly, the resultant magnetic interaction between these 
constituents is also changed appreciably. It is proved as the result of general 
study that the anomalous increase in thermo-remanent magnetization of ferro- 
magnetic minerals is due to magnetic coupling of a constituent having a higher 


Curie-point with other one having a lower Curie-point. 
1. Introduction 


It is an established fact that the thermo-remanent magnetization of ferro- 
magnetic minerals in rocks generally decreases monotonously with increasing temper- 
ature when heated in non-magnetic space [1], [2]. But some cases have been found 
where the T.R.M. increases anomalously around the temperature just below Curie-point 
during the process of thermal-demagnetization [3], [4]. This particular phenomenon in 
the temperature change of T.R.M. may be named the anomalous increase of T.R.M.. 
This phenomenon should be distinguished from the Hopkinson effect: the latter is 
the increase of the induced magnetization, instead of the remanent one, just below 
Curie-point when heated in a constant weak magnetic field. 

Previously, one of the present authors (T.N.) has suggested the two following 
possible mechanisms to explain this anomalous increase of T.R.M., [3]; z.e. 

(i) P-type of Ferrimagnetism. 

(ii) Negative magnetic interaction between two different constituents contained 

in the ferromagnetic minerals. 

In the first explanation, the anomalous increase of T.R.M. is attributed to a 


* Contribution from Division of Geomagnetism and Geoelectricity, Geophysical Institute Tokyo 
University. Series II. No. 63. ; 
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particular type of temperature change of spontaneous magnetization of ferrimagnetics 
named P-type by L. Néel [5]. 

On the other hand, the second explanation assumes the existence of negative 
magnetic interaction between two magnetically different constituents in the ferromagne- 
tic minerals, both having an ordinary type of the temperature change of spontaneous 
Magnetization, but with different Curie-points. The basic idea of this interpretation is 
quite similar to that applied for interpreting the phenomenon of the reverse thermo- 
+ remanent magnetism [6], [7], [8], [9]. The 

second mechanism is illustrated schematically 
in Fig 1, where the curves A and B show 
the temperature change of TR.M. of the 
A assumed two constituents, the Curie-points of 
which are 9, and 6, respectively. 

Now, if there is a negative interaction 
between these two constituents during the 
process of weak field cooling from a _ tem- 
perature higher than 9,4, the magnetic do- 
mains of the constituent having the lower 
Curie-point 9, may be, at its domain fixing 
temperature, affected by the demagnetizing 
Fig. 1 Schematic representation of anomalous field owing to the magnetization of a pre- 

increasing phenomenon of T.R.M. 
viously fixed constituent having a_ higher 


@, = % 


Curie-point 6 4. 

Consequently, if this demagnetizing field acting on the constituent with a lower 
Curie-point is strong enough to overcome the externally applied field, and to fix the 
domains of this constituent in the reverse direction, the resultant T.R.M. should be the 
difference between Curves A and B, which is shown by the dotted curve C in Fig. 1; 
the dotted curve C shows the anomalous increase of T.R.M. which we are going to 
deal with. 

In this paper, the mechanism of the above mentioned anomalous phenomenon of 
T.R.M. of the magnetite sand of Nii-sima derived from biotite rhyolite pumice will be 
mainly discussed. 

In the following discussion, we shall refer to the constituent with the higher 
Curie-point as H-constituent and that with the lower Curie-point as L-constituent. The 
result of the present studies seems to prove that the concerned anomalous increase of 
T.R.M. is due to the negative magnetic coupling between H- and ,L-constituents, as 


presumed in the previously proposed hypothesis (ii). 
2. Thermo-Magnetic Studies 


As stated in the introduction, it has previously been reported that the tempera- 
ture change of T.R.M. of magnetite sand of Nii-sima Island shows the anomalous 
increasing phenomenon during the process of its thermal demagnetization; the T.R.M. 
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decreases at first with increasing temperature, but increases abruptly near the Curie- 


point [3]. 

In order to study this particular phenomenon the magnetite sands were collected 
from Habuse-Ura in Nii-sima Island of Izu-sichito. 

First the main results of the present experiments will be summarized as follows. 
The original sample reveals no anomalous increase of T.R.M., but it begins to appear 
after a heat-treatment at a high temperature, say 550°C, in air. The amount of the 
anomalous increase is nearly proportional to the length of time during which the 
sample was kept at the high temperature. (Exp. A-(a)~A-(e)). 

The Curie-point of the original sample is about 470°C (Fig. 2), but after heat- 
treatment at 550°C in air for 23 hours, the sample comes to possess two Curie-points 
at about 470°C and 540°C (Fig. 5), indicating that the sample now contains two 
magnetically different constituents. These two constituents correspond to the Z- and 
H-constituents mentioned in the last paragraph of 1. 

On the other hand, when we heat the sample which has been sealed in a quartz 
tube evacuated to as low as 10-*mmHg, neither the anomalous increasing pheno- 
menon of T.R.M., nor the appearance of the higher Curie-point can be observed (Fig. 7). 

Furthermore, if the heat-treatment in air is made more thoroughly, raising the 
temperature of heat-treatment up to 650°C, the sample becomes a single constituent 
(or H-constituent), and at the same time, loses the characteristic of the anomalous 
increase of T.R.M. (Fig. 9). Next, we shall describe our experimental procedures and 
results in rather more detail. 

(i) Sample: 

The sample consists of the grains of titano-magnetite which have been magne- 
tically separated out from the magnetite sand, their mean diameter being about 400/. 
Its chemical composition is FeO : 52.24, Fe,O; : 39.52, TiO,: 8.24 in mol percent accord- 
ing to K. Katsura. 

(ii) Thermo-magnetic experiments : 

The experiments are composed of the measurement of (1) thermal change of 
magnetization (yH—T) in H=20e shown by full circles in Fig. (2)~Fig. (10), (2) the 
mode of development of T.R.M. in H=20e shown by hollow circles, and (3) its ther- 
mal decay in null magnetic space shown by broken lines. 

The measurements were made by the ballistic method, by using a similar 
apparatus to the one described in Fig. 2-2 (Rock Mag. [3)]). 

(ii-A) Production of anomalous increase of T.R.M. by heat-treatments at 550°C in air: 

(a) Fig. (2) shows the above mentioned characteristics of the original sample; 
the T.R.M. decreases monotonously and disappears at 470°C. In this case, no anoma- 
lous increasing phenomenon is observed. 

Next the original sample was heated up to 550°C, kept at this temperature for 
2, 10 and 23 hours and cooled rapidly out of the furnace. Fig. 3, 4 and 5 illustrate 
the xH—T and the T.R.M.-characteristics of these heat-treated samples. 
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Fig. 2 Magnetic characteristics of the origi- 


nal sample. 


Full circles : Temperature change of magneti- 


zation ; yH—T in H=2 Oe 


Hollow circles: Development of total T.R.M. 


in H=2 Oe 
Dotted line: Thermal decay of T.R.M. 
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Fig. 3 Magnetic characteristics of the sam- 

ple after heat-treatment at 550°C in air for 

1 hour. The T.R.M. remains nearly constant 

from 470°C to 500°C in the process of ther- 

mal decay. 

Full circles: Temperature change of magneti- 
zation; yH—T in H=2 Oe 


Hollow circles: Development of total T.R.M. 
in H=2 Oe 
Dotted line: Thermal decay of T.R.M. 


(b) After heat-treatment at 550°C in air for 1 hour: The T.R.M., as shown by 
the thermal decay curve of T.R.M. in Fig. 3, decreases monotonously up to 450°C, and 
remains nearly constant from 470°C to 500°C, above which it decreases again to dis- 


appear at 540°C. 

(c) After heat-treatment at 550°C in air for 10 hours: The T.R.M. decreases 
monotonously up to 450°C, and increases abruptly at this temperature as shown in 
Fig. 4. 

(d) After heat-treatment at 550°C in air for 23 hours: The thermal decay 
curve of T.R.M. in null magnetic field reveals, as shown in Fig. 5, a remarkable 
anomalous increasing phenomenon at 480°C. In this case, even in the process of 
development of T.R.M., a peak corresponding to the anomalous increasing phenomenon 
can be detected. 

In the yH—T curve in the heating process, two distinct Curie-points at about 470°C 
and 540°C can be identified. This indicates the existence of two constituents in this 
heat-treated sample. Judging from the experimental results described above, it may 
be reasonably considered that the constituent with Curie-point at 540°C, or the H- 
constituent, has been developed by the heat-treatment. 

(e) In Fig. 6, the change of partial thermo-remanent magnetization : J3n vy 20e 
with increasing temperature in null field for four samples which underwent the heat- 


treatment at 550°C in air for 0.5, 2, 10 and 23 hours respectively, and for the original 
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0 100 200 300 400 500 600°C 


Fig. 4 Magnetic characteristics of the sam- 
ple after heat-treatment at 550°C in air for 
10 hours. The T.R.M. decreases monoto- 
nously up to 450°C and increases abruptly 
at this temperature in the process of thermal 
decay. 


Full circles: Temperature change of magneti- 


zation; yH—T in H=20e 

Hollow circles: Development of total T.R.M. 
in H=2 Oe 

Dotted line: Thermal decay of T.R.M. 
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Fig. 5 Magnetic characteristics of the sample 
after heat-treatment at 550°C in air for 23 hours. 
Thermal decay curve shows a remarkable anoma- 
lous increasing phenomenon at 480°C, even in the 
process of development of total T.R.M., a peak 
corresponding to the anomalous increasing pheno- 
menon can be observed: in the yH—T curve, two 
distinct Curie-points at about 470°C and 540°C 
can be identified. : 
Full circles: Temperature change of magnetiza- 
tion; yH—T in H=2 Oe 
Hollow circles: Development of total T.R.M. in 
H=2 Oe 
Dotted line: Thermal decay of T.R.M. 


sample are shown. It will be clearly seen that the longer the time of heat-treatment, 
the more marked the anomalous increasing phenomenon, and that the temperature at 
which the anomalous increase attains its maximum value is kept nearly constant, 
namely about 490°C independent of the duration of heat-treatment. 

(f) In the above experiments: (a)-(e), all the heat-treatments were made in the 
atmosphere. In order to evaluate the role of oxidation during the heat-treatment, we 
made the similar heat-treatment in vacuum, and examined whether or not the same 


effect would also take place. 


In this heat-treatment, the sample was sealed in a 


quartz-tube evacuated to as low as 10-*mmHg. 
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Fig. 6 Thermal change of partial T.R.M. 

(J3s0-c, #=20e) With increasing temperature. 

Fig. 6-a: Original sample. 

Fig. 6-b: After heat-treated at 550°C for 
0.5 hour in air. 

Fig. 6-c: After heat-treated at 550°C for 
1 hour in air. 

Fig. 6-d: After heat-treated at 550°C for 
10 hours in air. 

Fig. 6-e: After heat-treated at 550°C for 
23 hours in air. 


Fig. 7 illustrates the results of measurement made on the sample heat-treated at 
550°C for 2 hours in the evacuated quartz-tube. Both the xH—T curve and the develop- 
ment of the total T.R.M. are similar to these of original sample shown in Fig. 2, and 
no anomalous increasing phenomenon is found in the thermal-demagnetizing curve of 
T.R.M. 

(ii-B) Effect of heat-treatment at 650°C in air: In this case, the temperature 
of heat-treatment was raised to 650°C, and the same measurements were made, #.e. ; 
(1) the temperature change of magnetization : xH—T in H=2 O8, (2) the development 
of total T.R.M. in H=2 Oe, and (3) the thermal demagnetization of the T.R.M. in null 
magnetic field. 

The results are similar to those of the experiments in (A) as long as the dura- 
tion of heat-treatment does not exceed some ten hours, but with longer heat-treat- 
ments, different features begin to appear in the thermal decay mode of the T.R.M. in 
null field; the anomalous increasing phenomenon disappears again. 
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0 100 200 300 400 = 500 600°C 0 100 200 300 400 500 600°C 
Fig. 7 Magnetic characteristics of the sam- Fig. 8 Magnetic characteristics of the sam- 
ple after heat-treatment at 550°C for 2 hours ple after heat-treatment at 650°C in air for 
in evacuated-quartz tube. Both yH—T and 3.5 hours. 
the mode of development of the total T.R.M. Full circles: Temperature change of mag- 
are similar to these of original sample (see netization ; yH—T in H=2 Oe 
Fig. 2) and no anomalous increasing pheno- Hollow circles: Development of total T.R.M. 
menon is found. in H=2 Oe 
Full circles: Temperature change of mag- Dotted line: Thermal decay of T.R.M. 


netization; yH—T in H=2 Oe 

Hollow circles : Development of total T.R.M. 
in H=2 Oe 

Dotted line: Thermal decay of T.R.M. 


(a) After heat-treatment at 650°C in air for 3.5 hours: (Fig. 8) 

In the xH—T curve, the Curie-point of both constituents, 470°C and 540°C are 
evident and the anomalous increasing phenomenon appears at about 480°C in the ther- 
mal decay curve of T.R.M. A slight peak can be detected in the curve of the develop- 
ment of total T.R.M. at about 450°C which may be considered as corresponding to the 
anomalous increase in the thermal decay curve of the T.R.M. 

(b) After heat-treatment at 650°C in air for 12 hours: (Fig. 9) 

The xH—T curve shows that this sample consists of the H-constituent nearly 
exclusively. The T.R.M. decreases monotonously with increasing temperature, and 
disappears at 540°C without any anomalous increasing phenomenon. This fact may 
be explained as follows: the time of heat-treatment is in this case long enough for the 
oxidation to invert the original constituent z.e., the L-constituent, into the H-constituent 
completely, so that there remains no L-constituent to interact with. 
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Fig. 9 Magnetic characteristics of the sam- 
ple after heat-treatment at 650°C in air for 
12 hours. The xH—T curve shows that this 
sample consists of the H-constituent nearly 
exclusively and the T.R.M. decreases mo- 
notonously with increasing temperature 
without any anomalous increasing pheno- 


menon. 


0 100 200 300 400 500 600°C 


Fig. 10 Magnetic characteristics of the sam- 
ple previously heated at 550°C in air for 10 
hours and then heat-treated at 1100°C for 
4 hours with subsequent quenching in a 
quartz-tube evacuated to as low as 10-*mmHg. 
Full circles: Temperature change of mag- 
netization; yH—T in H=20e 
Dotted line: Thermal decay of T.R.M. 


Full circles: Temperature change of mag- 
netization ; yH—T in H=2 Oe 

Hollow circles : Development of total T.R.M. 
in H=2 Oe 

Dotted line: Thermal decay of T.R.M. 


(ii-C) Other heat-treatment ; Disappearance of anomalous increasing phenomenon 


due to the homogenization by quenching: 
From the above experiments (A) and (B), it has been concluded that the samples 


having the characteristic of anomalous increase of T.R.M. are composed of two 


magnetically different constituents. 
To obtain another verification of this conclusion, we examined the effect of 


quenching from high temperature on the anomalous increasing phenomenon; if the 
temperature of quenching is high enough for this two constituent system to be homo- 
genized, the anomalous increasing phenomenon will no longer be observed. 
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For this purpose, heat-treatment at 1100’C for 4 hours and subsequent quench- 
ing were made for the sample sealed in the quartz-tube evacuated to as low as 
10-*mmHg. The sample sealed in the quartz-tube had previously been heat-treated at 
550°C for 10 hours in air to acquire the anomalous increasing characteristics. (See, 


Exp. A (c)) 
In Fig. (10), the xH—T curve, in H=2 Oe, and the thermal decay mode of T.R.M. 


of this quenched sample are shown. 

The results are as follows; the decay mode of T.R.M. shows no anomalous in- 
creasing phenomenon, and the xH—T curve indicates that the sample, now, has a single 
Curie-point at 510°C, which is just between the Curie-points of H- and L-constituents. 
Furthermore, it may be noticed that the intensity of the T.R.M. at room temperature 
is increased 4~5 times by this quenching (cf. Fig. 2-9). This increase of the intensity 
of T.R.M. at room temperature may also be explained by the homogenization of the 
two constituents: for there are no longer negatively fixed domains. 

An X-ray analysis of this quenched sample by Norelco spectrometer also proved 
that the sample consists of a single constituent with the crystal parameter; a=8.403 
+0.0018A, which is intermediate between those of constituent Z and H. (See section 
3 below) 


3. Chemical and X-ray Analyses 


In the previous sections, it has been pointed out that the anomalous increasing 
phenomenon has been attributed to the magnetic interaction between two different 
constituents Z and H, one of which has been developed by oxidation of the other 
constituent (Z). 

In order to examine the above-mentioned two-constituent model more precisely, 
chemical and X-ray analyses have been carried out. The results of these examinations 
obtained in each sample are summarized below. The chemical composition has been 
determined by T. Katsura of the Dept. of Geochemistry of Tokyo Institute of Tech- 
nology. 


Table I 
A > Se yee ) ingen ag Seager oh | . | 
Sample : _in Mol. Crystal | Custeonel | 
| -point 
| FeO | Fe, TiO, Parameter | eg 
(i) Original | 52.24 | 39.52 | 8.24 | 8.40840.0009A| 470°C 


(ii) Heat-treated for 23 hours 4g 79 | 44.736, 3.55 | 8-395+0.001A | 540°C 


in air at 550°C 8.404+0.002A 470°C 
(iii) Heat-treated for 12 hours 
poe sierted i 43.88 | 47.40 | 8.71 et 540°C 


In the X-ray analysis of sample (ii) as well as in the xH—T curve shown in 
Fig. 5, the existence of the two phases corresponding to the constituents Z and H is 
evident, while in the data for sample (i), which shows no anomalous increasing pheno- 
menon, only one phase is seen. 


Although sample (iii) showed no anomalous increasing phenomenon, and was seen 
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to be nearly a completely single 
phase H in the xH—T curve in 
Fig. 9, very faint peaks (a= 
8.404 A) corresponding to the 
constituent Z were observed, 
together with the strong peak of 
the H-constituent in the X-ray 
spectrograph. 

Therefore, we may expect 
that such a small amount of the 


constituent Z can not result in 
any noticeable anomalous increas- 


ing phenomenon. In all cases, 
no other phases than the above 


mentioned ones, were observed, 


; a. Fig. 11 Chemical composition of various heat-treated 
and chemical composition of the samples represented on a FeO-Fe,0,-TiO, diagram in 
moll %. 


constituent H may be considered 


“ . 5a 1: Original sample. 
as being quite similar to that of 2: After heat-treated at 550°C in air for 2 hours. 
sample (iii). In the ternary dia- 3: After heat-treated at 550°C in air for 10 hours. 
gram below (Fig. 11), are shown 4: After heat-treated at 550°C in air for 23 hours. 
i oA 5: After heat-treated at 650°C in air for 12 hours; this 
the chemical compositions of composition is nearly similar to that of the consti- 
either original and heat-treated tuent H. 


samples. 


4. Microscopic Examinations 


In this section, the results of microscopic examination of the polished surfaces 
of ferromagnetic mineral samples dealt with in the previous sections will be described. 

The surfaces of examined samples were all etched with HF for 5~6 seconds. 
By these examinations, a distinguishable correspondance between the occurrence of the 
anomalous increasing phenomenon and the development of a peculiar laminated struc- 
ture was verified. 

In the polished surface of the original sample which showed no anomalous 
increasing phenomenon, no intergrowth was observed (Fig. 12), but in those of the 
sample which had undergone the heat-treatment in air (heated in air at 550°C for 0.5, 
10, 23 hours) to acquire the anomalous increasing characteristics, a well developed 
laminated structure was observed. (Fig. 13-15) 

These lamellae were seen to have been etched more deeply than the host. Since 
it was proved that similar intergrowth of ferromagnetic minerals can originate the 
magnetic interaction that produces the Reverse T.R.M. [12], we may also reasonably 
assume that these lamellae will be the very structure producing anomalous increasing 


phenomenon, provided that the host and the lamellae consist of the Z-constituent and 


H-constituent respectively. 
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Fig. 12 Photo-micrograph of polished surface of Fig. 13 Photo-micrograph of polished surface of 
the original sample, which shows no anomalous the sample heat-treated at 550°C in air for 0.5 
increasing phenomenon, after etching with HF hour, after etching with HF for 6 seconds: a 
for 6 seconds: no intergrowth can be observed. faint laminated structure can be found near the 
grain boundary. 


Fig. 14 Photo-micrograph of polished surface Fig. 15 Photo-micrograph of polished surface of 

of the sample heat-treated at 550°C in air for the sample heat-treated at 550 C in air for 23 

10 hours, after etching with HF for 6 seconds: hours, after etching with HF for 6 seconds: the 

the well developed laminated structure can be laminated structures spread over the whole sur- 
observed especially at grain houndary. faces. 


Moreover, the following facts may be noticed in Fig. (13)-(15); (i) the laminated 
structure seen in the polished surfaces of these heat-treated samples is more remarkable 
as the time of heat-treatment in air increases, (ii) the said structures, begin to develop 
at the outer boundaries of the grains (Fig. (13), (14)), and spread over the whole sur- 
faces as the time of heat-treatment in air increases. (Fig. 15) 

For more precise examination of these laminated structure, the grains showing 
a marked anomalous increasing phenomenon were examined with the aid of an elect- 
ron-microscope, the observation being carried out by the triahol-alminum replica 
method. 

Fig. 16 is the electron-micrograph of the etched surface of the sample heat- 
treated in air at 550°C for 10 hours. The width of these lamellae is approximately 
about 0.017, while the length ranges between 0.5/¢ and St. 
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Fig. 16 Electron-micrograph of polished surface of the sample 
heat-treated at 550°C in air for 23 hours, after etching with HF 
for 5 seconds. 


5. Detailed Mechanism of the Anomalous Increasing Phenomenon of T.R.M. 


From the preceeding experiments, it has become clear that anomalous increasing 
phenomenon is originated by the negative magnetic interaction between two magnetic- 
ally different constituents, Z and H, the latter constituent being produced by oxidation 
to form the laminated structure. There have been several investigations reporting that 
when more than two magnetically different constituents are present in materials, and 
the configuration of their co-existence is in a laminated structure, negative magnetic 
interaction exists between them [12], [13]. 

In the following, the mechanism of the anomalous increasing phenomenon will 
be discussed in detail, and the results of the two experiments illustrating the above 
mentioned magnetic interaction, will be demonstrated. 

Now, let 7,” and T,% denote the domain fixing temperatures of the constituent 
H and L respectively with a reasonable assumption that 7,” > T),7.. When the sample 
is cooled in a magnetic field 4») from a temperature higher than 7,”%, the domains 
of the constituent H are fixed at 7,” in the direction of the applied field 2,,»p result- 
ing in an irreversible magnetization J/”(T) for T<T,”, then the effective field heg on 


the domains of the constituent Z at T= 7,’ bocomes, 
Neg =Ntapp—-N- JT). at T=T,* (1) 
where, N is the demagnetizing factor determined by the nature of the spatial configu- 


ration between the constituent Z and H. 
The necessary condition for the production of the anomalous increasing pheno- 


menon, then, will be just as in the case of the reverse T.R.M. [7], 
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(i) Field dependency of the anomalous increasing phenomenon appearing in the 
thermal decay of T.R.M. 

Fig. 17 shows the demagnetization of T.R.M. produced in various magnetic fields, 

by heating in a non-magnetic space. 


COEF 


In these figures, the manner in 
which the anomalous increasing pheno- 
menon ceases to appear as the intensity 
of the applied field is increased, will be 
observed. This disappearance of 
anomalous increasing phenomenon with 
the increase of the applied field may 
be interpreted as follows; when Japp is 
increased to an intensity of, say 2008, 
the right-hand side of Eq (1) becomes 
positive, so negative intraction exists 
no longer. 

(ii) Field dependency of the 

anomalous increasing pheno- 


menon appearing in thermal 


0 100-200 300 400 500 600°C or 
; ema .T.R.M. 
Fig’ 17 Thermal decay of T.R.M. (Jz,,,) produced . a = gnetization of P.T.R.M 
in various magnetic ficid ‘h). It can be seen that Use 
the anomalous increasing phenomenon ceases to In the last paragraph, we have 
a AS RES : ‘ ; : ; 
<r e intensity of the applied field (4) is jointed out that anomalous increasing 
eurve-Az Srcacines curve. BY” Jncidns phenomenon appearing in the thermal 
Curve. © * Sacissons CUrVes Ds Jnakcinne demagnetization of T.R.Ms. ceases to 


curve E: Jren-sooe Curve F: Jr, r-sove occur if the intensity of the applied 


field exceeds a critical value, say, 20 Oe. (Fig. 17) But the nature of the field depend- 
ency of anomalous increasing phenomenon appearing in the thermal-demagnetization. 
of partial T.R.Ms was found to be quite different from that of the total T.R.Ms men- 
tioned above. (5 (i)). 

First some theoretical considerations with regard to this point will be given 
here, and their results will be compared with those of the experiments. 

In Fig. 18, is shown the thermal-demagnetization curve of the P.T.R.M.; J rs yr Pro- 
duced by field cooling through the temperature range; 9, to Ty, and subsequent 
cooling in zero magnetic field, where 6, is the Curie-point of the constituent H and 
T, is to be taken as being a little higher than the Curie-point of the constituent ZL. 
(9,) Now, let us define 47 shown in Fig. 18, as the measure of the intensity of 
anomalous increase. 

The effective field eg on the magnetic domains of the constituent Z at T= yi 
will be readily seen by taking /,»»=0 in Eq (1), since the applied field is absent at 
this temperature. 


feg= N-J"(T,") (3) 
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On the other hand, since 4/, the inten- 
sity of anomalous increase, may be 
considered to be proportional to this 
effective field (3), we obtain, 


4I <x heg=—N-J"(T,*) 


ru or 4[ co J"(T,*) (4) 


The right hand side of Eq (4), namely, 
the irreversible magnetization of the 
G4 Gp constituent H; J"(T,4) may be calcu- 
Fig. 18 Definition of AJ as the measure of the lated as follows. 
intensity of anomalous increase. 

As has been stated already, the 
constituent H produced by oxidation of the constituent Z has a laminated structure in 
the matrix of the constituent Z. (Fig. 16) Now, for simplicity in calculation, we con- 
sider lamellae as being thin square plates whose dimensions are estimated from the 
electron-micrograph, (see Fig. 16) as follows; the width is approximately 0.01 and the 
side is of the order of lw. This size may be considered to be sufficiently small to 
give a single domain structure [14], and consequently, when the external field % is 
applied parallel to this plate, the change of magnetization proceeds only by the rota- 
tion of the total magnetic moment m of this plate. 

On the other hand, the volume of each of these lamellae being so small (v= 
10-'*cm, see later), it is readily supposed that the effect of thermal fluctuation with the 
order of magnitude of kT/vJ can not be negligible, when the temperature is raised at 
about 500°C. 

Therefore, the magnetization of the plate, parallel to # will be, using the Boltz- | 
man statistics, 


m=vJ,-tanh ( aaa a (5) | 
where J,#(T) stands for the spontaneous magnetization of the constituent H at 7. 
Since, at the temperature T= 7", the magnetic domains of the constituents 1 
are fixed, the irreversible magnetization of the constituent H of which the demagneti- 
zing field acts on the magnetic domains of the constituent L becomes approximately 
at T=T,4, (the increase in the spontaneous magnetization of H from Ty" ‘to 7T,* is 


: 
assumed negligible), : 
(Ti!) ) 
i) na 3 : 


In(T,) « tanh ( PTH 


and from Eq (4), 


| 
H H h 
AI « tanh (= eed ) (6) 


In actual case, the angle between the direction of # and that parallel to the ; 
plates is not fixed, so being averaged out for the whole orientation of the plates, 4/ | 


becomes, 
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pi ong. ("tanh an cos 0} -cosd do (7) 
7 vf b 


Now, we may assume that 7,”~=500°C, and J.(T,")=5, for the present sample, 
and that the volume of this lamellae as, v=(10~4)?x 10-6=107'*(cin)’, from the electron- 
micrograph, shown in Fig. 16. 

Inserting these values into the right-hand side, we obtain the curve of Eq. (7) vsh 
by graphical integration as illustrated. by the full line in Fig. 19, where the value 
of Eq. (7) is taken as 1 for h=oo. In 
this figure, it is remarkable that this 
fuuction (7) (or theoretically estimated 


4%/aJ,.. 


4I) approaches a constant value for 
the field intensity of several Oersteds. 

Now, let us compare these results 
with the experimental ones. In fact, 
such a saturation effect of 47 was also 
observed in the experiments. We 


examined the anomalous increasing Op 20 2D» © SD 26 ©. Bw 


phenomenon in the thermal demagneti- Fig. 19 Plot of the observed values of AJ,/AI,... 


zation of Partial T.R.Ms: J2°@,, pro- against the applied magnetic field & (hollow circles), 


; B and its theoretical curve calculated from Eq. (7). 
duced in various field /, for the sample 


heat-treated in air at 550°C for 10 hours. The observed value of 4J vs # are plotted in 
Fig. 19, where abscica denotes the ratio (47),/(4D),-soce- 

From this figure, it is obvious that 4/7 is saturated to a constant value when the 
intensity of the external magnetic field becomes of the order of magnitude of 10~20 Oe. 
It may be concluded that the theoretical results are in fair accordance with the experi- 
mental ones. 


6. Conclusive Remarks. 


When two or more different ferromagnetic minerals are contained in rocks, there 
generally exists some magnetic interaction between them. Several examples have been 
reported in which such magnetic interaction brings about certain anomalous pheno- 
mena of geophysical importance. 

For example, the Reverse T.R.M. found in the pumice of Haruna, the pitchstone 
of Asio, and the magnetite sand of Sokota has been found to be originated by magne- 
tic interaction between two ferromagnetic members belonging to the ilmenite-hematite 
solid solution [12], [15]. 

In this paper, we have investigated another example of the magnetic interaction 
between two members of the titano-magnetite solid solution series ; we have confirmed 
that the anomalous increasing phenomenon appearing in the thermal-demagnetization 
of T.R.M. is due to a kind of magnetic interaction between two constituents H and L. 

Through the magnetic and crystallographic studies, it was revealed that their 
Curie-points are 470°C and 540°C, and that both of them are of the spinel structure 
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with the lattice parameter a=8.404 A and a=8.395 A. 

It will be important to notice that the constituent H can be produced by oxida- 
tion and developed as lamellae in the matrix of the constituent L as. was observed 
by the microscope; in other words the formation of the constituent H is caused not 
by exsolution of a homogeneous phase of the solid-solution, but by oxidation of the 
constituent ZL. 

Finally, it will be also worthwhile to note that the point representing the cons- 
tituent H-spinel (cf. 3) lies midway between Titanomagnetite series and Ilmenite-Hema- 
tite series, since this part of the ternary diagram seems not to have been studied well 
so far. [16], [17], [18]. 

In conclusion, the writers wish to express their sincere thanks to Messrs. 
Akimoto and Uyeda, who helped them throughout the whole course of the present 
study. Their cordial thanks are also due to Prof. H. Kuno for his continued coopera- 


tion on the petrological side. 
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Horizontal Wind Systems in the Ionospheric E Region Deduced 
from the Dynamo Theory of the Geomagnetic Sz Variation 


Part I. Non-rotating Earth 


By Hiroshi MAEDA 
Geophysical Institute, Kyoto University 


(Read Oct. 29, 1955; Received May 30, 1956) 
Abstract 


Horizontal wind systems in the ionospheric E-region are obtained from the 
dynamo theory of the solar daily magnetic variation on quiet days, for the case 
of a non-rotating earth. It is found that the main harmonicterm of ionospheric 
wind velocities necessary to produce the quiet-day magnetic variations is not 
semi-diurnal, but diurnal. The diurnal amplitude is about 30m/sec and the 
semi-diurnal one is about 10 m/sec at low latitudes. The semi-diurnal component 
is nearly in phase with that of winds on the ground derived from the barometric 
pressure variation, and its tidal amplification is estimated as about 50. On the 
basis of these ionospheric winds, the vertical drift velocities of electrons in the 
F-region are estimated. 


1. Introduction 


In recent years several radio and non-radio methods have been developed [1] for 
measuring the velocities of horizontal movements taking place in the ionosphere, but 
there is a question as to whether the movements detected by the various radio tech- 
niques are true movements of the air. 

According to the dynamo theory, put forward by Stewart [2], the geomagnetic 
variations are due to electric currents in the upper atmosphere, produced by motion 
of electrically conducting air across the geomagnetic field. Recent results obtained in 
the United States [3] by means of rockets suggest strongly that these electric currents 
flow mainly in the E-region of the ionosphere. Hence it will be possible to deduce 
ionospheric winds from the geomagnetic variations by the dynamo theory, if the dis- 
tribution of ionospheric conductivity is known. A treatment along this line was first 
made by Vestine [4], who deduced upper-atmospheric winds from the dynamo theory 
of geomagnetic disturbance. 

Very recently, the present author [5] has proposed a formula for the world dis- 
tribution of ionospheric conductivity on quiet days. By using this formula, the wind 
systems in the ionospheric E-region required from the dynamo theory of the quiet-day 
daily magnetic variation (S,) are obtained. Part I deals with a simple case in which 
it is supposed that the motion of the air is horizontal and irrotational. 


(12h) 
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2. Method of Treatment 


According to the dynamo theory the electric field E is the sum of the electro- 
static field H, and the “dynamo” electric field E;, that is, 


E.=E,,+Eiz; (4) 
E,=EytEiy ’ 


where the subscripts x and y refer to the southward and eastward components, res- 
pectively. Let S be the electrostatic potential, then 


os 
i; = ’ 
aoe 
(2) 
ie hee 
sis asin 002’ 


where a is the radius of a current sheet (about 6.48x108cm), @ the colatitude, and A 
the east longitude. Since we supposed that the motion of the air is irrotational, the 
velocities are derived from a potential @ as follows: 


“= — ne 
aoe’ 
3 
ole aed es 
asin 00,’ 
and the dynamo fields are given by 
E;, ren vH, ’ ( 4 ) 
Evy pal me uH,, 


where H, is the vertical component of the geomagnetic field. Assuming for simplicity : 
that the geomagnetic field is that of the centred dipole, alined along the geographical | 


axis, we can write 


H,=—Ccos@, 


where C is a constant (about 4 gauss Hence, from (1), (2), (3) and (4) we have 


7 fe eects, 
00 Oa (5) 
0S 00 
et mere nie Cod 0 
<n ta | geal”) 
Eliminating S from (5) gives 
ed ek as 00 yatke ae (E sind} 
cotd + 50 (sin dcosa5) = C aa: z) ag °" 
=F(0, 4). (6) 


This is the differential equation to be satisfied by the velocity potential @. We deal | 
with here only the periodic components of mean annual winds, hence the following . 


boundary conditions must be taken: 
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@=0 at the poles (for @=0 and @~=z), 
7 
a =0 at the equator (for a= >) ba 


3. Solution of the Differential Equation 


A solution of the differential equation (6) may be found in the following manner. 
If we assume that @ is expressed as a series of surface harmonics, that is, 


O= 331 (k™ cos mA+ Rnysin mA) P (cos 9), (8) 


nm 


where P™(cos@) denotes the Schmidt’s normalized spherical functions [6], and 4 the 
local time reckoned in angular measure, then the left-hand side of (6) can be trans- 


formed [7] into the form 


sin 0 3133 (k%,cos mA+ ke ,sinmd) G2 (8), (9) 
where 
G0) = — ene. {{(0-1)?— 1] [e+ 1)2— me}? Pm, , + [22 em}? Ps). 


The right-hand side of (6) is calculated from the distribution of the total electric field 
E. This field has recently been obtained by Maeda [5] from the relation 


1( ~ 5, )=BieEt Deny ’ 
2n 


1,( = =a = —DiskK.+ DwEy ’ 
2n 


using the geomagnetic north (4X) and east (4Y) components of S, corrected by C, 
(because of this correction, the zonal term of EH almost vanishes) and the height-inte- 
grated conductivities given by the form 


per OS Gos DSvy = Dive ® (X), Day = Die? (2) 


where ¥(x)=1+2.00cos x + 1.46 cos* x, 
and x=solar zenith angle. 
The results are shown in Figs. 1A and1B. The right-hand side of (6) calculated by 
using this electric field may also be expressed as a series of surface harmonics, that 
is, 
F(6, 4)=sin 6 f (8, 2) 
=sin 93331 (f¥,acos MA +f X sin Ma) P¥ (cos 8), . (10) 


where the values of coefficients /%’s are shown in Table 1. By equating the corres- 
ponding terms of (9) and (10), we have 
Rr aGn(O) =f¥ Px (cos 0), 
Kit ,G3(0) = f% PX (cos0). a 
It is expected from the coefficients of P%,, and P™_, in (9) and the values of 
¥’s in (10) that the harmonic coefficients k%’s for each order decrease fairly rapidly 
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Table 1. Spherical harmonic coefficients of 


| >| 
SIN | 
PR £ F(8, 2)= a> (f¥,qc0s M+ f¥ ,sin M)) PM(cos 6). 
2 ———- 8 (Unit : 10!°cm?/sec) 
= — g -—— — ve ———— 
Load | y } 
am Serta i. 1 2 327.54 244.67 
an se y 4 152.05 42.16 
tA D 6 47.01 185.19 
em é 8 81.88 77 82 
A - — ~— —_ = aes 
- 2 | s —176.21 —92.80 
00 04 ' 08 12 16 20 py ‘ot —76.51 49.47 
LOCAL TIME (HouRS) 7 117.57 105.03 
Fig. 1A. Distribution of the total electric field. 9 182.85 124.71 
3 4 94.71 62.49 
6 | —166.60 78.43 
8 —48.27 46.21 | 
10 53.96 53.04 


Fig. 1B. Latitude-distributions of the harmc 

nic amplitudes A’s and phases 0’s (or f,,’s, th 

time at which the maximum of the mth terr 

occurs) for the total electric field component 
E,, Ey. 
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with increasing degree. We can therefore evaluate approximately the coefficients ky’s 
from (11) by neglecting higher terms. The result of calculation is shown in Table 2. 
And with these values the expression (8) gives an approximate solution of the differ. 
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Table 2. Spherical harmonic coefficients of 
® (0, 2) = (k™ cos m).+ k™ , sin m?) P™(cos 9). 
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(Unit: 10'°cm*/sec) 
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Fig. 2B. Distributions of_the wind velocity for each of the 24-, 12- and 8-hourly components. 
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ential equation (6). This solution is a particular solution of (6), and it is easily 
seen that this particular solution satisfies the boundary conditions (7: The differ- 
ential equation (6) is linear. Its general solution is given by the sum of a particular 
solution and the complementary function. As will be seen in Appendix, however, the 
complementary function of (6) does not satisfy the boundary conditions. It is therefore 
concluded that the particular solution obtained here is the unique solution of (6) to 
satisfy the boundary conditions (7). 


4. Wind Systems 


By using the spherical harmonic coefficients of ® shown in Table 2, the velocity 
is calculated from relation (3). The latitude distributions of the harmonic coefficients 
for the southward and eastward velocity components are shown in Fig. 2A, and the 
distributions of winds are illustrated in Fig. 2B for each of the 24-, 12- and 8-hourly 
components. The total wind system of these components is shown in Fig. 2C. As will 
be seen in these figures, the main harmonic term of wind velocities is diurnal. The 
diurnal wind has a source (at about 15”) and a sink (at about 3”) on the equator, and 
its amplitude is about 30m/sec. The semi-diurnal wind has two sources (at about 1.5” 
and 13.5") and two sinks (at about 7.5" and 19.5"), and its amplitude is about 10 m/sec 
at low latitudes. 
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Fig. 2C. Distribution of the total wind velocity. 
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It may be of interest to compare these results with the wind velocities on the 
ground derived from the pressure variation. Among the 24-, 12-, 8- and 6-hourly 
components of the solar barometric variation, the regularity of the 12-hourly (semi- 
diurnal) component is most remarkable. Hence we here take up only the semi-diurnal 
component. This component was considered, in detail, by Simpson [8]. If we take 
the semi-diurnal variation as given by the first term of Simpson's. formula 

p=0,937 sin (2 wt+154°) mmHg, 
where ¢ is local time and » the angular velocity of the earth, then the southward and 
eastward compoents of the air velocity on the ground, derived from this pressure 
variation, are obtained [9] as follows: For a non-rotating earth 


ut» = 32.8sin*0 cos @ cos (2 wt+154°)cm/sce, v,= —21.8sin20sin (2 wt 4-154") cm/sec, 
and for the rotating earth 
2 = 54.6 cos 8 cos (2 wt+154°)cm/sec, v,= —(21.8+32.8 cos*@) sin (2 wt + 154°) cm/sec. 


At the equator, for both these cases, #,=0 and the amplitude in v, is about 22cm/sec. 
Comparing these winds with the ionospheric winds obtained above, we find that the 
semi-diurnal component of ionospheric winds is nearly in phase with that of winds on 
the ground, and that its tidal amplification is about 50. 


5. Vertical Electron Drifts in 
the F-Region 


Using the ionospheric winds ob- 
tained above, we can easily calculate 
the dynamo field from (4). The results 
are illustrated in Figs. 3A and 3B. 
And, by subtracting the dynamo field 
from the total electric field, we have the 
electrostatic field as shown in Figs. 4A 
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Fig. 3B. Distribution of the dynamo field. 


Fig. 3A. Latitude-distributions of the harmo- 
nic amplitudes A’s and phases 0’s (or t»’s) for 
the dynamo field components E;,, Eiy. 
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N LATITUDE Ss 
Fig. 4A. Latitude-distributions of the harmonic 
amplitudes A’s and phases 6’s (or ¢,’s) for the 
electrostatic field components E,,, Egy. 
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Fig. 5. Daily variations of the vertical drift 
velocity of electrons in the ionospheric F-region 
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Fig. 4B. Distribution of the electrostatic field. 


and 4B. 


electrons in the F-region from the relation 


case in which no tidal oscillation exists in the 
F-region (full lines), (2) the case in which the 
same tidal oscillation exists in the F-region as in 
the E-region (dotted lines), and (3) the case in 
which the opposite tidal oscillation exists in the 
F-region to in the E-region (dashed lines). 
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W= Fr cos 


[10], where ¢ is the magnetic dip angle, and F the intensity of the geomagnetic field 
in the F-region. Since there is some uncertainty about winds in the F-region, calcula- 
tions are made for the following three cases: (1) the case in which no tidal oscillation 
exists in the F-region, (2) the case in which the same tidal oscillation exists in the 
F-region as in the E-region, and (3) the case in which the opposite tidal oscillation 
exists in the F-region to in the E-region, and the results are shown in Fig. 5. 

As will be seen in the figures, the diurnal component of the drift velocities is 
predominant. The time of maximum upward drift is about 8" at latitudes higher than 
20°, but it becomes about 11’ at the equator. This phase (11’?) of the vertical drift 
at the equator is very convenient [11] to interpret the observed daily variations of 
maximum electron density N,, and its height #,, in the equatorial F2-region at sun- 
spot-minimum. For the phase (8"t) at latitudes about 20° the effect of the drift on 
N,,F 2 is very slight [12] and moreover the ion production in these latitudes is roughly 
the same as at the equator. It is therefore possible to interpret the fact that N,,F2 
at about 20° is greater than that at the equator. Thus the main feature of the geo- 
magnetic distortion in the F2-region near the equator (between latitudes about + 20°) 
can be interpreted by the effect of the vertical electron drifts. 


6. Concluding Remarks 


The wind systems which cause to give rise the geomagnetic S,-variations are 
obtained from observed geomagnetic data by the dynamo theory of anisotropic con- 
ductivity of upper layers. As original data are the total S,-field, including the part 
due to currents induced in the earth’s interior, the results obtained here are overesti- 
mated just by this internal part. In order to compare the present results with observ- 
ed ionospheric winds, it must also be remembered that the data refer to the mean 
solstice, %(S+W), say. In recent years, several methods have been developed for 
measuring the ionospheric winds. Although the results of these direct measurements 
are not always regarded as to give true movements of the air, it may be noted that 
the wind velocities obtained by our calculation are the same order of Magnitude as 
those observed. 

The electric field EZ was estimated by Maeda (see [5]), by eliminating the effect 
of the daily variation of conductivity from S,, under the assumption that the S,-elect- 
ric currents flow in one layer (the E-layer). It was found that the main harmonic 
term in E is diurnal. Consequently, also in the wind velocities obtained here, the 
diurnal component is most predominant. The diurnal amplitude is about 30 m/sec, and 
the semi-diurnal one is about 10m/sec at low latitudes. If we compare this result 
with winds on the ground derived from the barometric pressure variation, it is found 
that the semi-diurnal component of ionospheric winds is nearly in phase with that of 
winds on the ground, and that the tidal amplification is about 50. This ratio should, 
however, be reduced by the reason mentioned above. Fejer [13] calculated the elect- 


130 H. Marepa 


ric current-densities in the ionosphere from the assumption of wind system by the 
dynamo theory, and estimated a tidal amplification of about 60. His result, however, 
is clearly an overestimate, because he compared the electric currents produced only by 
the semi-diurnal winds with the total S,-currents. 

The semi-diurnal component of ionospheric winds will be explained by the 
theory of atomospheric oscillations [14]. The origin of the diurnal component is not 
yet clear, but it may be suggested that this is produced by any thermal process. 

On the basis of the ionospheric winds obtained here, the vertical drift velocities 
of electrons in the ionospheric F2-region were calculated, and this result supports a 
view [11, 12] that the main feature of the geomagnetic distortion in the F2-region 
near the equator can be interpreted by the effect of the vertical electron drifts accom- 
panied by the S,-currents. 

This paper dealt with a simple case where the motion of the air is irrotational. 
A further study in which the hydrodynamic equations on the rotating earth are taken 
into account will be made in Part II by Kato [15]. 
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APPENDIX 


Complementary Function of the Differential Equation for 0 
In order to obtain the complementary function of the differential equation (6) for 


@, we must solve the homogeneous equation 
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If we put 
0(6, 4)=X(9)Y(A), 


then (A-1) can be separated into the following two equations: 


lhe aX) _ nrcot0 X=0 A-2 
fj (sind coso<, m*cot : (A-2) 
ay. mY=0, (A-3) 


dt 
where m? is a constant of separation of variables. The linearly independent solutions 
of (A-3) are easily found as follows: 


Y,=cos mA, Y,=sin mA. 


To obtain the solution of (A-2), put 
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then (A-2) becomes 


z(1—z) oe + [(m+1)—(m+2)z] = (+ 1 =o. (A-4) 


This is the hypergeometric equation; the cases of m=1 and m=3 being non-degene- 
rate cases, and the case of m=2 being a degenerate case, and its solution is well 
known [16]. Thus two lineary independent solutions of (A-2) about =0 are obtained 


as follows: 


: | | 
m First solution | Second solution | 
| OEE YS See SR iF FO FES 
1 X,= sin oF(,. 9 pagahae) | X, = X,log(sin*6) + sin OF (,. 9 + 2.5 sin 0) | 
2 X,=sin?0F(1, 2 ; 3; sin?9) X,=1/sin*9 
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We must examine whether or not these solutions satisfy the boundary conditions 
(7), The second solutions X, tend to infinity as 0-0, so that these solutions do not 
satisfy the condition at the pole. On the other hand, the values of dX,/d@ tend to 
infinity as 0-7/2, so that the first solutions do not satisfy the condition at the equator. 


We can therefore conclude that the complementary function of (6) does not satisfy the 
prescribed boundary conditions. 
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